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Abstract 

Quantized Skyrmions with baryon numbers B = 1,2 and 4 are considered 
and angularly localized wavefunctions for them are found. By combining 
a few low angular momentum states, one can construct a quantum state 
whose spatial density is close to that of the classical Skyrmion, and has 
the same symmetries. For the B = 1 case we find the best localized 
wavefunction among linear combinations of j = ^ and J = | angular 
momentum states. For B = 2, we find that the j = 1 ground state 
has toroidal symmetry and a somewhat reduced localization compared 
to the classical solution. For B = A, where the classical Skyrmion has 
cubic symmetry, we construct cubically symmetric quantum states by 
combining the j = ground state with the lowest rotationally excited 
J = 4 state. We use the rational map approximation to compare the 
classical and quantum baryon densities in the B = 2 and B = 4 cases. 
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1 Introduction 



The connection between the quantum and classical descriptions of a many- 
body system is an important but rather tricky one. In nuclei, the existence 
of a rotational band, a sequence of states whose energy increases with angular 
momentum j approximately as + 1), where / represents a moment of 

inertia, suggests the existence of a static intrinsic classical shape to the nucleus 
which is not spherically symmetric [1 . It is not obvious how this classical shape 
arises, and it is hard to predict the shape, but one can partially reconstruct it 
from the spectrum. 

For a rigid body, the quantum states of various angular momenta are given 
by Wigner functions Di^{a, (3, 7), where a, P, 7 are the Euler angles parametriz- 
ing the orientation, j is the total angular momentum and s, m its components 
with respect to the body-fixed and space-fixed third axis. Symmetries of the 
body constrain the possible s-values or combinations of s-values that can oc- 
cur. A classically oriented state is a (5-function in the Euler angles. This can 
be obtained by taking an infinite linear combination of Wigner functions. For 
a body with symmetry, one would take a sum of 5-functions on a set of orien- 
tations related by symmetry (which are not distinguishable). Even if there is 
no fundamental rigid body to start with, one can consider these linear combi- 
nations. Thus, given a rotational band of states, one can construct a classically 
oriented state by taking an infinite linear combination of true quantum states 
of definite angular momentum. The properties of this oriented state (e.g. the 
particle density) would define the nature of the intrinsic state. 

Something like this has been done in certain condensed matter situations. 
One may construct a classically oriented state when all that is rigorously avail- 
able is quantum states labelled by angular momentum. Cooper et al. j2] have 
studied a model of rotating states of a Bose condensate trapped in a harmonic 
well (whose shape essentially makes the condensate two-dimensional). By nu- 
merically combining precise states over a range of angular momenta, they have 
shown that a condensate with vortices can be obtained. In the rotating frame, 
these vortices form a static array, and so are angularly localized despite the 
rotational invariance of the problem. The vortices do not really exist in any of 
the states of definite angular momentum, but they do in the combined state, 
and they can also be physically observed. 
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This localization depends on the system being large. Ideally, the moment of 
inertia should be almost infinite. In that case, the angular momentum states of 
different j are almost degenerate, and the angular localization may be achieved 
at almost no energy cost. (Similarly, an object with large mass can be spatially 
localized by taking a superposition of momentum states.) 

Unfortunately, for nuclei, this is not always a realistic way to proceed. For 
large nuclei, like Hf^^°, there are many states in a rotational band, and it is 
pretty clear that an intrinsic nuclear shape exists P]. For smaller nuclei, how- 
ever, at most a few low-lying states can be identified as forming a rotational 
band, and their energy separation is quite large because the moments of iner- 
tia are smaller. Not much is known about the wavefunctions of the states in 
the band, so it is hard to consider linear combinations. Instead it is better to 
postulate some intrinsic shape and fit its parameters to data. In this way, it 
is found, for example, that the Ne^° nucleus has a prolate deformation but 
one cannot say it is exactly a prolate ellipsoid. 

An alternative treatment of many-body systems can give angularly (and 
spatially) localized states much more easily. This is the approach based on 
an effective field theory, for example a Ginsburg-Landau description of a Bose 
condensate. Here, classical solutions of the field equation can naturally exhibit 
spatial order, for example an array of vortices. Because of the underlying sym- 
metries, the classical solution is not unique, but is parametrized by collective 
coordinates describing, say, the center of mass position and angular orientation. 
Comparison with the previous discussion suggests that effective field theory can 
only be valid for large systems of many particles. To reconstruct quantum states 
of definite angular momentum, one may quantize the collective coordinates; this 
makes sense if the mass and the moment of inertia are of finite, but not infi- 
nite magnitude. A critical comparison of exact quantum states and classical 
solutions of an effective field theory has been carried out for quantum Hall 
ferromagnets by Abolfath et al. |3]. 

In this paper, we shall consider the Skyrme model and its connection with 
nuclei and their various angular momentum states. The Skyrme model is an 
effective field theory of pions, with a topological quantum number that can be 
identified with baryon number 0. Skyrme's original idea was that the model is 
justified because nuclei can be thought of as made up of a condensate of many 
light pions (with a topological winding). Recently, the justification is based on 
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the idea that each nucleon is made of Nc quarks, where SU{Nc) is the gauge 
group of QCD, so a nucleus of baryon number B is made of a large number, 
NcB, of quarks [Hj- The Skyrme model has a semi-rigorous standing if Nc is 
large, but it is a controversial matter whether the physical value A^c = 3 is 
sufficiently large. 

The classical Skyrme field equation, like that of the Ginsburg-Landau model, 
can be solved numerically and much is known about its minimal energy solutions 
(especially for pion mass equal to zero) 0. Most importantly, the classical 
shapes of the solutions, and their symmetries, are known for values of B up to 
and beyond 20 (and work is underway to take account of the finite pion mass, 
which could have a qualitatively significant effect for B > 10). These classical 
shapes could represent the intrinsic shapes of nuclei of modest size. 

The shapes obtained have no obvious relation to shapes of nuclei as under- 
stood using other models, in particular, models based on point nucleons. For 
example, four-nucleon potential models are used to describe the a-particle, and 
the classical minimum occurs for a tetrahedral configuration of the nucleons [S]. 
In the Skyrme model, the solution of minimal energy with B = 4 has cubic 
symmetry. 

Our aim in this paper is to bridge the gap between the classical Skyrmion 
shapes and the quantum states of nuclei. The traditional approach has been to 
quantize the collective coordinates of Skyrmions, seek the lowest energy states 
consistent with the allowed values of the angular momentum, and compare 
with the ground state properties of nuclei. This approach has some success in 
reproducing the known spins of nuclei, especially for even baryon number. More 
recently, a table of allowed angular momenta for the ground and first excited 
states of rotationally quantized Skyrmions has been constructed [H]. However, 
in these quantum states of definite angular momentum, the original Skyrmion 
shape information is sometimes completely lost. 

We cannot consider an infinite linear combination of angular momentum 
states, as we expect that large angular momenta will lead to Skyrmion de- 
formations, or if these are suppressed, then to infinite energy. Instead, here 
we shall consider an intermediate picture. By taking a small combination of 
low-lying angular momentum states, we partially reconstruct the shape of the 
classical Skyrmion solution. We shall optimize the angular localization of the 
Skyrmion within the limited combinations of states at our disposal. Such a 
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finite combination of states has finite energy (not necessarily very much higher 
than the ground state). If one could investigate theoretically (or experimentally, 
although this could be difficult) the same combination of angular momentum 
states in another nuclear model, one might see better the connection with the 
Skyrme model picture. 

For the B = 1 case we apparently do not have the problem of orientation 
because a single Skyrmion has a spherically symmetric density. However, the 
Skyrmion still has rotational collective coordinates, and we will show that a 
particular combination of j = | and J = | states gives the most localized wave 
function. We also show that the ground state of the deuteron (the j = I state 
of the B = 2 Skyrmion), without an admixture of higher angular momentum 
states, retains the toroidal symmetry of the classical solution. Forest et al. have 
argued that not only the pure deuteron state, but also deuteron clusters within 
larger nuclei, show toroidal structure jTHj. 

Finally, we shall show that a combination of j = and j = 4 collective 
states of the B = 4 cubic Skyrmion gives a state close to the classically oriented 
Skyrmion. The same combination of j = and j = 4 states in a four-nucleon 
potential model could be compared. Now it is inevitable that in the potential 
model, the state will have cubic symmetry because of the Wigner functions 
involved. In that sense our discussion is purely kinematic. However, some 
detailed properties of the state (density, currents) might show a close similarity 
with the Skyrme picture. 

This paper is restricted to the B = 1, 2, 4 cases, and is organized as follows. 
Section 2 contains a review of the Skyrme model (for more details see ^2)- 
Section 3 we give an outline of the rational map approximation for Skyrmions, 
and its consequences, and recall the rational maps for the B = 1, B = 2 and 
B = 4 Skyrmions. In Section 4, we use the j = | and j = | quantum states 
of a i? = 1 Skyrmion introduced in ^21, and find the combination which gives 
the best localized wave function. In Sections 5 and 6 we use the rational maps 
introduced in Section 3 to find the "best" wavefunctions for i? = 2 and B = A 
Skyrmions, respectively, and calculate the quantum baryon density in these 
states. In Section 7 we briefly discuss the implications of adding vibrational 
modes, and summarize our conclusions. 
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2 The Skyrme model 



The Skyrme model j3] is an effective low energy theory of QCD attempting to 
treat pions, nucleons and nuclei. The topological soliton solutions arising from 
this model can be interpreted as baryons. 
The model is defined by the Lagrangian 

^Trid.Ud^U^) (1) 



where U{t, x) is an S'f/(2)-valued scalar field. and e are parameters which 
can be scaled away by using energy and length units of -F^/4e and 
respectively. Thus, with the values of F^, and e as in [12] our units are related 
to conventional units via 

F 2 

— = 5.58 MeV , = 0.755 fm . 

4e eF^ 

Introducing the su(2)-valued right current i?^ = {d^U)U'^ and using geo- 
metrical units, the Lagrangian (1) may be rewritten in the concise form 



\Ti{R,Rn + ^Tr([i?^,i?,][/?^i?'^])| d'x. 



(2) 



The Euler-Lagrange equation which follows from Q is the Skyrme equation 

d,(R^ + \[R,,[R''.R']]^=^- (3) 

Static solutions are the stationary points (either minima or saddle points) of 
the energy function 



^ = Y^j {-^Tr(i?ii?,)-^Tr([i?„i?,][i?„i?,])|c?3x^ 



(4) 



where we have introduced the additional factor l/127r^ for convenience. 

U, at fixed time, is a map from M.^ into S^, the group manifold of SU{2). 
However, the boundary condition U 1 implies a one-point compactification 
of space, so that topologically U: ^ S^, where the domain is identified 
with M'^ U {oo}. As the homotopy group 773(6''^) is Z, maps between 3-spheres 
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are indexed by an integer, which is denoted by B. This integer is also the degree 
of the map U and has the exphcit representation 



As -B is a topological invariant, it is conserved under continuous deformations 
of the field, including time evolution. This conserved topological charge Skyrme 
identified with baryon number. 

Static fields of minimal energy, solving the Skyrme equation, are called 
multi-Skyrmions (Skyrmions, for short). They have been constructed numer- 
ically for B up to 22 7J, and the symmetries of these solutions have been 
identified. For B = 1 the Skyrmion has spherical symmetry, and for i? = 2 
toroidal symmetry. It turns out that Skyrmions have non-trivial discrete sym- 
metries for B > 2. Solutions for negative B are obtained by the transformation 
U —>■ W, which preserves the energy. 

3 Rational map ansatz 

In what follows, we will be using the rational map approximation to Skyrmions 
|14j . Rational maps were first introduced into the theory of three-dimensional 
solitons by Jarvis p3|, in the context of monopoles, but they prove to be very 
useful for Skyrmions as well. 

Rational maps are maps from S"^ —>■ S"^, whereas Skyrmions are maps from 
— S''^. The idea in [T^ is to identify the domain S"^ of the rational map with 
concentric spheres in M^, and the target of the rational map S"^ with spheres of 
latitude on S*^. A point in can be parametrized by (r, z); r denotes radial 
distance and the complex variable z specifies the direction. Via stereographic 
projection z can be written in terms of usual polar coordinates 6 and (p as 
z = tan(^/2)e*'^. A rational map may be written as R{z) = p{z)/q{z), where 
p{z) and q{z) are polynomials in z. The degree of the rational map, A^, is 
the greater of the algebraic degrees of the polynomials p and q. N is also the 
topological degree of the map (its homotopy class) as a map from S'^ ^ S"^. 
The point z on S"^ corresponds to the unit vector 




Similarly, the value of the rational map R is associated with the unit vector 

1 



Hi? 



;{R + R,i{R- R),l- 



(7) 



1 + 

The ansatz for the Skyrme field, depending on a rational map R{z) and a radial 
profile function /(r), is 

U (r, z) = exp{if{r)hR^^) ■ r) , (8) 
where r = (ri,r2,r3) denotes the triplet of Pauli matrices, and /(r) satisfies 

/(O) = TT, /(OO) = 0. 

An SU (2) Mobius transformation of z corresponds to a rotation in physical 
space; an SU{2) Mobius transformation of R (i.e. on the target 5*^) corre- 
sponds to an isospin rotation. Both these transformations of a rational map are 
symmetries of the Skyrme model, and both preserve A^. 

It can be verified that the baryon number for the ansatz (jSJ is given by 



B = - 



f f sin / 



27r2 



1 + 


z 


2 


dR 


1 + 


R 


2 


dz 



2i dzdz 
Jl + Wf 



r dr . 



(9) 



2idzdz/{l + is equivalent to the usual 2-sphere area element sm6d6d(j). 

The angular part of the integrand, 

(10) 



1 + 


z 


2 


dR 


1 + 


R 


2 


dz 



1 + k 



2^2 



is precisely the pull-back of the area form 2idRdR/{l + |-Rp)^ on the target 
sphere of the rational map R, so its integral is An times the degree N of the 
map. Therefore (jnj simplifies to 

-2N 



B 



TT 



f sin^ fdr = N. 



(11) 



An attractive feature of the rational map ansatz is that it leads to a simple 
energy expression which can be separately minimized with respect to the ra- 
tional map R and the profile function / to obtain close approximations to the 
numerical, exact Skyrmion solutions, and having the correct symmetries. (The 
numerical solutions are in fact best found by starting from the optimal rational 
map approximations.) 

Indeed, using (jH)) we get the following expression for the energy (jlj): 

-sin"^/' 



E 



1 
3^ 



r'f" + 2N sin' f if" + 1)+ 1- 



dr . 



(12) 
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Here X denotes the purely angular integral 



1 = ^ 

An 



1 + 


z 


2 


dR 


1 + 


R 


2 


dz 



2idzdz 

(13) 



1 + b 



2^2 



which only depends on the rational map R. To minimize E, for maps of given 
degree N, one should first minimize X over all maps of degree N. Then the pro- 
file function /, minimizing the energy (fT^ . may be found by numerically solving 
a second order, ordinary differential equation with and X as parameters. 

For B = 1, the rational map is R{z) = z, and this reproduces Skyrme's 
hedgehog ansatz 0, which is exactly satisfied by the B = 1 Skyrmion. For 
B = 2 and B = A the symmetries of the computed Skyrmions are -Doo/i and 
Oh respectively, and in each case there is a unique rational map of the desired 
degree with the given symmetry, which also minimizes X. They are, respectively, 

R{z)=z\ Ri^z) = ^—j=—— . (14) 

In all these cases, we have made a convenient choice of orientation in presenting 
the maps. 

When quantizing the Skyrme field, we will be interested in the behavior of 
the wavefunction with respect to different orientations of the Skyrmion con- 
figurations. Consequently, all the information we need will be encoded in the 
angular dependence of the baryon density ()10|). which only depends on the 
rational map, and the profile function / will not be of much interest for our 
purposes. 



4 B = 1 case 

The B = 1 Skyrmion is spherically symmetric and takes the hedgehog form 

f/o(x)=exp{2/(r)x-r)}, (15) 

where /(O) = n and f{oo) = 0. If Uq is the soliton solution, then U = AUqA~^, 
where A is an arbitrary constant SU (2) matrix, is a static solution as well. But, 
in order to get solitons which are eigenstates of spin and isospin one needs to 
treat A ClS db collective coordinate. So substitute 

U = A{t)UoA-\t) 
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in the Lagrangian (1), where A{t) is an arbitrary time-dependent SU{2) matrix. 
The Lagrangian for A is |2] 

(16) 



L = -M + \Ti{doAdoA~^] 



where M is the sohton mass and A is an inertia constant which may be found 
numerically. 

The SU{2) matrix A can be written as A = ao + ia. ■ t, with + |ap = 1. 
In terms of = 0, 1, 2, 3) the Lagrangian (fT^ becomes 



L = -M + 2A^(a^) 



(17) 



and after the usual quantization procedure one gets the Hamiltonian 



8A ^ 

5=0 



g2 



:i8) 



Because of the constraint + |ap = 1, the operator Yll=o (•9^/<9a|) is to be 
interpreted as the laplacian on the 3-sphere. The wavefunctions can be 
expressed as traceless, symmetric, homogeneous polynomials in the a^. 
Using the isospin and spin operators 
1 



Jk 



2* I «o 



d 


d 


d 


dak 


ak 

dao 


^klm.O'l „ 

dam 


d 


d 


d 




dak 


dam 



(19) 



Adkins, Nappi and Witten have found the normalized wavefunctions for neu- 
tron, proton and A-resonances [IS^. The wavefunctions we require, having 
opposite I3 and J3 eigenvalues, are 



|A",s, 



TT 



TT 
7T 



[ao - msj , 



(ao + ias) 



V2 



(ao + ia3)(l -3{al + al)), 



TT 

(ao - ia3){l - 3{al + aj)) 

TT 



V2, 



TT 



[ao - tas) 



(20) 
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But none of these "pure" j = 1/2 and j = 3/2 states is the best locahzed 
wavefunction. This will instead be given by a superposition of the above states. 

If we could take into account an infinite number of angular momentum states 
the most localized wavefunction would be the Dirac delta function, which may 
be expressed in the following form 

<5(/x) = J2 (2J + , J = 0,1,1,1,..., (21) 

j 

where X"' (/^) is the character of the representation of dimension (2j + l), and Oq = 
cos /i. However, this wavefunction does not respect the Finkelstein-Rubinstein 
(FR) constraints 16^, which in the case of one Skyrmion requires that the 
wavefunction is antisymmetric under A —>■ —A, thus ensuring that the quantized 
Skyrmion is a fermion. The sum in (PT|) must therefore be restricted to half- 
integer values of j, giving the total |(5(/i) — 6{fi — vr)). 

For the SU (2) group, the representation matrices are matrices of Wigner 
functions, i.e. for each j 



/ D^-- ... D^. ■ \ 

n 3-3 



\ -33 -3-3 / 

In what follows we will be interested in the j = 1/2 and j = ?)/2 cases, for which 
the Wigner functions take a concise form in terms of Oq, . . . , a^. The character 
is the trace of the above matrix, consequently let us write down the diagonal 
elements: 



nV2 _ I • 

1/2,1/2 ~ '^0 + ^C^S ; 

nl/2 _ 

-'^-1/2,-1/2 ~ '^0 ~ '^'^3 , 

n3/2 _ ( 1 • n3 

-'^3/2,3/2 ~ 1^0 + f'O'-ij , 

-^1/2.1/2 ^ (ao + m3)(l - 3(ai + 02)) 

^^-112-112 = (cto - ^a3)(l - 3(ai + 02)) 



^%'/2.-3/2 = («o-^a3)'. (22) 



-1/2,-1/2 
)3/2 

-3/2,-3/2 

If we truncate the sum (PT|) at j = 3/2 we get the following candidate for a well 
localized (normalized) wavefunction, 

^(ao, ai, a2, 03) = (al - "^a^ . (23) 
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In terms of nucleon and A-resonance states this can be written as 

^ (|A-) - |A°) - |A+) + |A++) - - b))) , (24) 

with spins as in (f^UI) . A more general wavefunction of this type is 

-1/2 



^2 / 5 \ ' 

^(ao, ai, 02, as) = — f — + K + (afj + koo) . 



(25) 



The maximum magnitude of \1/ at Oq = ±1 occurs when k = —3/8, confirming 
that this is the best locahzed wavefunction. 

Another measure of how well the wavefunction is localized around Oq = ±1 
is given by the integral 

( — -i- K. -i- ] 

7r2 



l^^ + K + K^^ al\al + Kaol'^dn. (26) 



Here oq = cos /i and dQ = An sin^ /i dfi is the measure of integration. After an 
easy calculation, we find that this integral is maximal when k, = —1/4, which is 
close to the result we got before. One more wavefunction worth considering is 

which is as well localized as the one with k = —3/8 according to criterion (j2ni), 
and rather simpler. It is the following combination of nucleon and A states: 

^ (|A-) - |A°) - |A+) + |A++) -2V2z{\n)- \p))) . (27) 

These localized states are not physically important for isolated nucleons; 
however, they could be useful for modelling nucleons in interaction. Recent 
developments have shown that, for example, the deuteron is not formed from 
a proton and neutron only, but probably also contains some amount of A- 
resonances f7|nHl- Therefore, considering a superposition of states with differ- 
ent angular momenta is definitely physically meaningful. In ^H] the deuteron 
was modelled by a bound state of Skyrmions in the attractive channel, where 
the relative orientation of the Skyrmions was chosen to maximize the attraction 
at short range. Such states could be approximated by the combined j = 1/2 
and j = 3/2 states we have discussed here. The dependence of the force be- 
tween two Skyrmions on their relative orientation is the classical analogue of 
the tensor force between nucleons, and it appears to automatically lead to an 
admixture of a A-resonance component to each nucleon. 
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5 B = 2 case 



The B = 2 Skyrmion has Dooh symmetry and a toroidal shape |20[l2H l^. and is 
used to describe the deuteron. We take the symmetry axis to be the third body- 
fixed axis, and the Skyrmion to be in its standard orientation if this coincides 
with the third Cartesian axis in space. In the rigid body approximation to 
quantization, the wavefunction is a function only of the rotational and isospin 
collective coordinates. (We ignore the translational collective coordinates, and 
set the momentum to zero.) To make an appropriate quantization we have to 
impose FR constraints, which tell us that the ground state has the quantum 
numbers = (0,1), where i is the total isospin and j is the total spin. 

The wavefunction describing this deuteron state was obtained in [23j. Since 
i = 0, there is no dependence on the isospin collective coordinates, and the 
(normalized) state is 



Here a, /? and 7 are the rotational Euler angles, D^^^a, (3,'j) is a Wigner 
function, and m is the third component of the space-fixed spin. 

In the analysis was extended to include one vibrational mode of the 
system, allowing the B = 2 toroidal Skyrmion to separate into two B = 1 
Skyrmions. The wavefunction therefore includes a factor u{p), the radial part 
of the deuteron wavefunction, which satisfies a radial Schrodinger equation on 
the interval [pq,oo), where po corresponds to the toroidal configuration. Here, 
however, we consider only the rigid body rotational states, and their angular 
dependence. 

Since we are particularly interested in the spatial orientation, we treat states 
differing in m as different. The state we are looking for has to have the same 
symmetry properties as the classical solution. Consequently, the desired wave- 
function is 



which is axially symmetric both on the left and on the right (i.e. with respect to 
the body-fixed symmetry axis, and the Xs-axis in space). \l/ has its maximum 
magnitude at /3 = and /3 = vr, corresponding to the Skyrmion in its standard 
orientation, and turned up-side down, which is classically indistinguishable after 
an isospin rotation. 






(29) 
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Now, given the orientational quantum state ()29|1 we may calculate the nu- 
cleon density, and find how quantum effects change the density of the classical 
configuration. We find the expression for the baryon density distribution (x) 
in physical space (which is interpreted as nucleon density) by averaging the 
classical baryon density over orientations weighted with I^E'p. 

The density in the quantum state is therefore 

p^(x) = J B{D{A)-^^)\'^{A)\^ sin (3 da dp d-f. (30) 

Here A stands for the SU (2) matrix parametrized by Euler angles a, f3, 7, and 
D{A) for the 5*0(3) matrix associated to A via 

D{A),, = ^Tiir^AnA^) . (31) 

As was already mentioned, the B = 2 rational map is R{z) = z'^, and this gives 
a good approximation to the B = 2 Skyrmion solution. It leads, using Q, to 
the following expression for the classical baryon density: 

1 /I + bp \^ 
i3(r,.) = -^^^|.|j ,(r), (32) 

where g{r) is a radial function. g{r) is unaffected by the quantum averaging, 
so we ignore it from now on. In terms of polar angles, the angular dependence 
of B is given by 



^ = Z^^^^^^TIRW^ ^ (33) 



1 (1 + tan2 (|)) tan^ 
^ (l + tan4(f))2 



where this is normalized to have angular integral equal to 2, the degree of the 
rational map. 

To evaluate p*(x) we first expand B in terms of spherical harmonics Yim{0, (p)'- 

B = Y,ClmYlm{eA), (34) 
l,m 

where, because of axial symmetry, there are only terms with m = 0. Although 
(|34j) is an infinite series it is a good approximation to take just the first two 
non-zero terms of the sum, 

i3 = CooFoo(^^) + C2oK,o(^), (35) 
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as all the other terms contribute less than a 5% correction. Because the map 
R has degree 2, cqo = l/v^; also, we find numerically that C20 = —0.36. Then 
B{D{A)^^^) can be written as 

B{±) = Cooroo(^~) + cA{e) , (36) 

where x = D{A)^^ x and similarly for 6, 0. Using the transformation properties 
of spherical harmonics under rotations, 

YimiO, 4>) = Y1 DlMyyikiO, 0) , (no sum on I) , (37) 

k 

the fact that = {3/8it'^)Dqq{A)DIq{A)* , the orthogonality properties of the 
Wigner functions 

/Q 2 
Di,iA)D'2Ar sin P da dp dj = ^-^d'^'SacSbd , (38) 

and (in terms of the Wigner 3j symbols) 

Di,{A)D{,{A)D{;{A)smPdadpdj = 8n'(^^^ ^'^ ^'^^ (^l 

•(39) 

we find directly from ()30p that the quantum probability distribution is 

2 

pq, = cqoYqo + -C20Y20 . (40) 
5 

This is an exact expression - no higher terms are present. We see that it 
resembles the classical distribution (j35|) . but is more dominated by the first 
term. Thus, when quantum effects are included, the classical toroidal density 
remains, but is smoothed out to become more spherically symmetric. 



6 B = 4 case 

The minimal energy 3 = 4 solution has cubic symmetry; the region of high 
baryon density resembles a rounded cube with holes in the faces and at the 
centre j21]. We define the orthogonal body-fixed axes to be those passing 
through the face centres, and the standard orientation of the cube to be where 
these axes are aligned with the Cartesian axes in space. We shall again consider 
the Skyrmion as a rigid body, which means the configuration is not allowed to 
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vibrate. It was shown in that in this case the ground state, representing 
the a-particle, has quantum numbers i = and j = 0, with the (unnormahzed) 
wavefunction \1>(°) = 1 being independent of the rotational and isospin collective 
coordinates. The first excited state has i = and j = 4: and is 



= Dt^ia, (3, 7) + ^J^D^c^, (3, 7) + Dt 



.4^(«,/?,7). (41) 



In |2n| the third component of the space-fixed spin, m, was arbitrary. The 
structure of (PT|l is required by the cubic symmetry with respect to body-fixed 
axes. 

But just fixing m is not enough to make the wavefunction cubically sym- 
metric both on the left and on the right, i.e. also with respect to space-fixed 
axes. To achieve this we need to take the following linear combination of the 
above wavefunctions: 

^(4) = vl/f + . + vl/W . (42) 
V 5 

The cubic symmetry in space is fairly obvious by analogy with (@H), and can 
be verified as follows. First note that symmetry under 90° rotations about the 
x^-axis implies that all possible terms in (j^^ with m other than ±4, vanish. 
To simplify the calculations a bit further we then introduce new variables 

a = cos(^^)e^^^e^^", 6 = - sin(^)e^^*^e^^" . (43) 

Obviously they satisfy |ap + |6p = 1. In terms of a and b, the SU{2) orientation 
matrix parametrized by Euler angles a,/3,7 is 



A 



a b 
—b a 



In this notation the wavefunctions for different m take the following compact 
form: 

= a« + 14aV + 6« 



^i'^ = (^aV + aV + ^(3 - 30(|a2| - I&'l)' + 35(|a2| - l^'l)')) 

= a^ + Ua%^ + b\ (44) 
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Therefore the wavefunction ()42|1 in terms of a and h is 



= 2Re(a« + 14aV + 6^) + 14(a^6^ + aV^ 

+ ^(3-30(|aVl&'lf + 35(|aVl&'ir) • (45) 
As expected, it is real. By acting on A with the generators of the cubic group: 



^j' l-¥ ¥ 



(46) 



corresponding to a 90° rotation around a face of the cube, and a 120° rotation 
around a diagonal of the cube, we find the resulting transformations of (a, 6), 
and it is easy to check that \E'^^-' is cubically symmetric both on the left and on 
the right. 

The wavefunction '^^'^^ has a positive maximum of y at the identity, (a, h) = 
(1,0), and at all other elements of the (double cover of the) cubic group. This 
is as desired, as it corresponds to the Skyrmion having a high probability to 
be in its standard orientation. But \E''^^^ also has a negative minimum of — 
which gives a further local maximum of |\E^'^^)p, at an orientation obtained by 
a 60° rotation around a diagonal of the cube, which is far from the standard 
orientation. We wish to suppress this. 

We can do this by being a bit more sophisticated than in the B = 2 case. 
We still have the freedom of adding an arbitrary constant to the wavefunction. 
This means taking a superposition of the ground and first excited states, \E'^°'' 
and 

^ = ^(4) + ft:^(o) . (47) 

Here again we are interested in the nucleon density of the configuration. Our 
goal will be to adjust the constant k to get a quantum distribution as close 
as possible to the classical one. As in the B = 2 case we define the quantum 
nuclear density via 



pvp(x) = J B{D{A)- 'x.)]"^ {A)\' sin f3 da dp d-f, (48) 

where B{x.) is the classical baryon density of the B = 4 Skyrmion in its stan- 
dard orientation. Using the rational map (fT^ . we find that B has the angular 
dependence 

B = -|z|2(l + \z\y ^ ^ ^ • (49) 
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Expressed in terms of polar angles, 



(tan«(f) -2tan4(f)cos40+ 1) 

X 



,2 ' 



(tan^d) + 2tan4(|)cos40 + 12tan4(|) + l)' 
which may be expanded in the following form: 

B = doYoo + d,z,{e, 0) + deZeie, 0) + dsZsie, </>) + . . . . (si) 

Here Z4, Zq and Zg are the unique cubically symmetric combinations of spherical 
harmonics with, respectively / = 4, 6 and 8:^ 

[U 

Z4 — Y44 + Y ~5~^4o + Y4-4 , 

[2 

— Yq4 — y -Yqq + Ye_4 , 



Zs — Yss + \ 777^84 + \ -^^80 + \ 777^8-4 + ^8-8 • (52) 
V 00 V DO V DO 

The leading coefficient is do = 2/1/7? because the rational map has degree 
4, and by numerical calculation we find that ^4 = —0.28, d^ = —0.032 and 
ds = 0.024. Then, by a similar calculation as in the B = 2 case, normalizing the 
wavefunction and using the orthogonality properties of the Wigner functions, we 
find the following numerical result for the angular dependence of the quantum 
baryon density: 

= c^o^o + TTTT^ ^{-(0.038 + 0.075/€)Z4 - O.OOGZg + 0.002^8) , (53) 

2.00 + I J 

which is again a finite sum, all the further terms being zero. 

In (jnni); 1^ is not yet specified. Let us adjust it in such a way that the above 
distribution looks as close as possible to the classical one, i.e. let us maximize 
the coefficient of the / = 4 terms: 

4 



2.5Q + K 



,.2 



(0.038 + 0.075k) 



^ These can be derived by combining the generating, cubically symmetric Cartesian poly- 
nomials a;^ + + , a;^ + y'* + , + + , and finding the combinations which satisfy 
Laplace's equation |27| . 
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The maximum is at k = 1.17, which leads to the following expression for the 
quantum baryon density: 

^ 1.13^00-0.13^4 - 0.006^6 + 0.002^8 

^ ^0^00 + 0.46^4^4 + 0.24^6 + 0.14^8. (54) 

Thus in the B — A case, as in the B — 2 case, one can find a quantum 
state which locahzes the Skyrmion close to its standard orientation, and which 
preserves the symmetry of the classical solution. However, the inclusion of 
quantum effects smoothes the classical baryon density, making it rather closer 
to spherically symmetric. Again, the effect is to approximately halve the leading 
non-constant harmonics, here with I — A. If we considered the pure j = 4 state 
we would get 

p^(4) ^ doYoo + 0.2^4^4 + 0.34^6 + 0.154^8 , (55) 

which is much closer to spherically symmetric. 
We can also find the energy of our state ^; it is 

E = ^ (2.56E,=4 + n^Ej=o) = 0.65E,=4 + 0.35E,=o , (56) 

so it is not as highly excited as a pure j — A state. 

The combination of j = and j — A states, ^, is a bit artificial as the 
quantum state of a free B — A Skyrmion, but would make sense if we were 
dealing with interacting Skyrmions (for example, when describing compound 
nuclei such as Be^, C^^ in the Skyrme model equivalent of the ct-particle model). 
Here we expect the relative orientations of the B — A subclusters to be rather 
precisely fixed when they are close together, so as to minimize their potential 
energy. 

7 Conclusions 

Three well-locahzed wavef unctions of the B = 1 Skyrmion have been considered 
and some of their advantages and physical implications have been discussed. 
The B — 2 and B — A minimal energy Skyrmion solutions have been quantized 
in such a way that the wavefunctions have the same symmetry properties as 
the classical Skyrmions (respectively, axial and cubic symmetry both on the left 
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and on the right), and angularly localized quantum states with shapes closest 
to the classical solutions have been found. In the B = A superposition 
of two low-lying states of definite angular momentum needed to be considered. 
It is impossible to completely reproduce the shape of the classical solution this 
way. The quantum state necessarily smoothes out the classical baryon density, 
making it closer to being spherically symmetric. 

All our results were obtained in the rigid body approximation, i.e. we did 
not allow the Skyrmions to vibrate. Considering the vibrational modes may 
be an interesting topic for future work. The vibrational modes for the B = 2 
and B = 4 Skyrmions were calculated in [2HI 121] and a qualitative analysis has 
been given for B = 7 [30]. The vibration frequencies obtained can be separated 
into those below and those above the breather mode, which is the oscillation 
corresponding to a change in scale size of the Skyrmion. In the modes below 
the breather, the Skyrmion tends to splits up into individual B = 1 Skyrmions 
or small B Skyrmion clusters. 

But treating the vibration modes as harmonic oscillators is not very accurate, 
since, as the minimal energy configuration separates into individual Skyrmions 
the potential flattens out. A more accurate treatment would involve estimating 
the inter-Skyrmion potential at intermediate and large separations. Thus it 
should not be expected that the inclusion of the zero point energy of harmonic 
vibrational modes will yield accurate results for masses, binding energies of 
states, etc. 

The vibrational modes also couple to the rotational degrees of freedom, 
which complicates the analysis of the rotational and isospin wavefunctions [31] . 
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